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The most recent results from direct searches for dark matter particles in the galactic halo are 
examined in terms of an effective Minimal Supersymmetric extension of the Standard Model at the 
electroweak scale without gaugino masses unification. We show that the annual modulation effect 
at 8.2 a C.L. recently presented by the DAMA Collaboration, as the result of a combined analysis 
of the DAMA/Nal and the DAMA/LIBRA experiments for a total exposure of 0.82 ton yr, fits 
remarkably well with what expected for relic neutralinos for a wide variety of WIMP distribution 
functions. Bounds derivable from other measurements of direct searches for dark matter particles 
are analyzed. We stress the role played by the uncertainties affecting the neutralino-quark cou- 
plings arising from the involved hadronic quantities. We also examine how present data on cosmic 
antiprotons can help in constraining the neutralino configurations selected by the DAMA effect, 
in connection with the values of the astrophysical parameters. Perspectives for measurement of 
antideuterons possibly produced in the galactic halo by self-annihilation of neutralinos belonging 
to the DAMA configurations are examined. Finally, we discuss how findings at LHC would impact 
on these issues. 
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I. INTRODUCTION 

We are at present witnessing a great activity in running 
and preparing experimental projects of direct searches for 
dark matter (DM) particles in the galactic halo [l[ . 

A number of new results came out recently. Most no- 
tably, the DAMA Collaboration, by analyzing the data 
obtained with the DAMA/LIBRA experiment Q, has 
confirmed the evidence for the annual modulation effect 
already measured with the previous DAMA/Nal exper- 
iment As a result of the combined analysis of the 
DAMA/Nal and the DAMA/LIBRA experiments, for a 
total exposure of 0.82 ton yr, the effect of the annual 
modulation is now at 8.2 a C.L. Q. Other Collabora- 
tions have reported upper bounds on the WIMP-nucleon 
scalar cross section 111, H i, Q (as usual here WIMP 
stands for a generic Weakly Interacting Massive Parti- 
cle), using other approaches and scaling from different 
nuclei. 

In this paper we discuss the relevance of these exper- 
imental results for the most widely discussed candidate 
for cold dark matter, the neutralino. Neutralino configu- 
rations (concerning here mainly light neutralinos) which 
fit the annual modulation data are extracted from the 
DAMA results and confronted with the current measure- 
ments on galactic antiprotons. Perspectives of seeing sig- 
nals due to these configurations in forthcoming measure- 



ments of galactic antideuterons are examined. The cur- 
rent experimental upper bounds on the WIMP-nucleon 
scalar cross section are also scrutinized. 

In all our considerations we explicitly keep into ac- 
count the uncertainties affecting the neutralino-quark 
couplings because of the involved hadronic quantities, 
typically the pion-nucleon sigma term a^N, the quan- 
tity Co which is related to the SU(3) symmetry breaking 
and a quark mass ratio r (see definitions later on). 

In the analysis of the DAMA data, we discuss both the 
case in which the channeling effect |7| is included in their 
analysis as the one where this effect is not taken into ac- 
count; these two cases are treated separately. We recall 
that the channeling effect occurs when an ion, travers- 
ing a crystalline structure along a path (quasi-) parallel 
to crystallographic axes or planes, has a deeper penetra- 
tion than in an amorphous material of the same compo- 
sition [8!]. In Ref. it is shown that the occurrence of 
the channeling effect makes the response of the DAMA 
Nal(Tl) detector more sensitive to WIMP-nucleus inter- 
actions than previously estimated. 

Finally, in the present paper it is stressed that CERN 
LHC (Large Hadron Collider), which hopefully will shed 
light on the Higgs sector and supersymmetric extensions 
of the Standard Model (SM), has a remarkable discovery 
potential in terms of the population of light neutralinos 
which fit the annual modulation data. 
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The plan of the paper is the following. In Sect. II the 
main features of the supersymmctric model employed in 
the present paper are delineated. Sect. Ill is devoted 
to a recollection of the main phenomenological proper- 
ties used in our analysis, and Sects. IV-V to a detailed 
comparison of our theoretical evaluations to the current 
results of WIMP direct detection. In Sect. VI we discuss 
the constraints implied by measurements of galactic an- 
tiprotons and in Sect. VII we examine the perspectives 
for forthcoming measurements on cosmic antideuterons. 
Impact of LHC on these issues is discussed in Sect. VIII, 
and finally conclusions arc drawn in Sect. IX. 



II. THE SUPERSYMMETRIC MODEL 



we take BR(B° S — > /j~+/j + ) < 1.2-10 -7 ; measurements of 
the muon anomalous magnetic moment a M = (g^ — 2)/2: 
for the deviation Aa^ of the experimental world aver- 
age from the theoretical evaluation within the Standard 
Model we use here the range —98 < Aa M ■ 10 11 < 565, 
derived from the latest experimental [l^ and theoretical 



17jj data. 



Also included is the cosmological constraint that the 
neutralino relic abundance does not exceed the maxi- 
mal allowed value for cold dark matter, i.e. £l x h 2 < 
(ficM/^lmax. We set (O c dm h 2 ) max = 0.122 as de- 
rived at a 2a level from the results of Ref. [18j. We 
recall that this cosmological upper bound implies on the 
neutralino mass the lower limit m x > 



7 GeV 0. 



The supersymmctric scheme we employ in the present 
paper is the one described in Ref. an ef- 

fective MSSM scheme (effMSSM) at the electroweak 
scale, with the following independent parameters: 
M\ , M2, p, tan /3, tua, Wq, mj and A. Notations are as fol- 
lows: Mi and M 2 are the U(l) and SU(2) gaugino masses 
(these parameters are taken here to be positive), p is the 
Higgs mixing mass parameter, tan (3 the ratio of the two 
Higgs v.e.v.'s, m,A the mass of the CP-odd neutral Higgs 
boson, m,q is a squark soft-mass common to all squarks, 
mr is a slcpton soft-mass common to all sleptons, and 
A is a common dimensionless trilinear parameter for the 
third family, A^ = Af = Arriq and A? = Amr (the tri- 
linear parameters for the other families being set equal 
to zero). In our model, no gaugino-mass unification at 
a Grand Unified (GUT) scale is assumed. The lightest 
neutralino is required to be the lightest supersymmctric 
particle and stable (because of R-parity conservation). 

In the present paper the numerical analyses are per- 
formed by a scanning of the supersymmctric param- 
eter space, with the following ranges of the MSSM 
parameters: 1 < tan/3 < 50, 100 GeV < \p\ < 
1000GeV,5GcV < Mi < 500 GeV, 100 GeV < M 2 < 
1000 GeV, 100 GeV < m^mj < 3000 GeV, 90GcV < 
m A < 1000 GeV, -3 < A < 3. 

The following experimental constraints arc imposed: 
accelerators data on supersymmetric and Higgs boson 
searches (CERN e+e" collider LEP2_[l(j and Collider 
Detectors DO and CDF at Fermila b II ill ): measurements 
of the b ~ ► s + 7 decay process [12J: 2.89 < B(b — > 
s+7)-10~ 4 < 4.21 is employed here: this interval is larger 
by 25% with respect to the experimental determination 
[l2| in order to take into account theoretical uncertainties 
in the SUSY contributions [l3[ to the branching ratio of 
the process (for the Standard Model calculation, we em- 
ploy the recent NNLO results from Ref. [Hj]); the u ppe r 
bound on the branching ratio BR(B® — > /i~ + /i + ) [l5j: 



III. PHENOMENOLOGY RELATED TO THE 
DIRECT SEARCH OF GALACTIC WIMPS 

In case of WIMPs with coherent interactions with nu- 
clei the differential event rate dR/dEn {En being the 
nuclear recoil energy) measured in WIMP direct searches 
can be written as 

dR pa rriTv . 2t . (nucleon) j-,2 1 771 \n-i \ n\ 

dE~ R = :Nt 7^2% A F (E R )J(v min ), (1) 

where: 

2(fmin) = / <TW . (2) 

Jw>Vmin W 

In these formulae, notations arc: Nj> is the number of 
target nuclei per unit mass, mjv is the nucleus mass, /tti 
is the WIMP-nucleon reduced mass, A the nuclear mass 
number, ct^^ 0011 ^ is the WIMP-nuclcon coherent cross 
section, F(Ejt) is the nuclear form factor, £ is the fraction 
of the mass density of the WIMP in terms of the total 
local density for non-baryonic dark matter po (i.e.: £ = 
Pw/po), /es(^) and w denote the velocity distribution 
function (DF) and WIMP velocity in the Earth frame, 
respectively (w = \w\), v m - ln is the minimal Earth-frame 
WIMP velocity which contributes to a given recoil energy 
E R : « min - [m N E R /(2^ A )} 1 / 2 ( MA being the WIMP 
nucleus reduced mass). 

It is convenient to define a velocity distribution func- 
tion in the Galactic rest frame f(v), where v = w + 
U0 being the Earth velocity in the Galactic rest frame. 
The Earth frame velocity DF is then obtained by means 
of the transformation: /es(^) = f(w + v^). By def- 
inition, the velocity distribution function f(v) is given 
by the six-dimensional phase-space distribution function 
F(r, v) evaluated at the Earth location R in the Galaxy, 
i.e. f(v) — F(Rq,v). The velocity DF f(v) is truncated 
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at a maximal escape velocity v osc , since the gravitational 
field of the Galaxy cannot bound arbitrarily fast WIMPs. 

By employing Eqs. ([T]) and ([2|) one can derive informa- 
tion on the quantity ^ct^^. 0011 ' from the measurements 
on the differential rate dR/dEn, once a specific velocity 
distribution function f(v) is selected. 



A. WIMP distribution functions 

A systematic study of various distribution functions 
was carried out in Ref. |l9| and a number of these DFs 
were subsequently analyzed in Ref. [2oT ] to study the sen- 
sitivity of the upper bounds on C^scaiar 00 "^ derived from 



different experiments of WIMP direct detection. 

In Ref. [l9| the phase-space DFs were classified into 
four categories, depending on the symmetry properties 
of the matter density (or the corresponding gravitational 
potential) and of the velocity dependence: A) spheri- 
cally symmetric matter density pdm with isotropic veloc- 
ity dispersion, B) spherically symmetric matter density 
with non-isotropic velocity dispersion, C) axisymmetric 
models, D) triaxial models 2l|, |22j. For each category, 
different specific models are identified. 

In Ref. a procedure was developed to determine, 
for each individual velocity DF, the relevant range of 
the pa values. Specifically, each halo model was con- 
strained by a number of observational inputs: i) proper- 
ties of the galactic rotational curve, namely the range 
of the allowed values for the local rotational velocity, 
170 km sec" 1 < v Q < 270 km sec -1 0, 0], and the 
amount of flatness of the rotational curve at large dis- 
tances from the galactic center, and ii) the maximal 
amount of non-halo components in the Galaxy, Af v j s 
(i.e. the disk, the bulge, etc.). These constraints de- 
termine the extremes of the local dark matter density 
po- For instance, when one assumes a maximal halo 
(i.e. the contribution of the non-halo components is min- 
imized) po is maximal and, on the contrary, the value 
of po is minimal when the contribution of the halo to 
the rotational curve is minimal. By this procedure, for 
any specific analytic form of the velocity DF and any 
value of the local rotational velocity vo within the inter- 
val 170 km sec -1 < vo < 270 km sec -1 , one derives the 
relevant lower and upper bounds, p™ 1 " and Po" ax , for the 
local DM density. Obviously, p™ in and Po" ax are increas- 
ing functions of «o, since a large amount of local matter is 
necessary to support a large value of the local rotational 
velocity. 

Notice that both p and Vq are crucial parameters in es- 
tablishing the properties of the quantity ^ct^^ , when 
this is extracted from the detection rate. Whereas po has 
the role of a normalization factor, vo is both related to 



the WIMP kinetic energy and to the change in the refer- 
ence frames, thus is crucial in determining the shapes of 
the sensitivity levels usually reported in the plane WIMP 

>. (nucleoli) 

In the present paper, as our reference model we take 
the cored isothermal sphere, which we will discuss in de- 
tail; we will also comment about some other DFs. The 
density profile of the cored-isothermal sphere (denoted 
as Evans logarithmic model, or Al model, in Ref. flit ]) 
is: 



p(r) 



3R,, 



4ttG (R 2 + r 2 )< 



(3) 



where G is the Newton's constant, vo is the local value of 
the rotational velocity and R c is the core radius. For R c 
we use the value R c — 5 kpc. For the parameter vo we will 
consider the values vo = 170,220,270 km sec -1 , which 
represent the minimal, central and maximal values of vo 
in its physical range. For these values of vo, the extreme 
values of po are: i) vq = 170 km sec -1 with p™ m = 0.20 
GeV cm -3 and p^ nax = 0.42 GeV cm -3 ; ii) v = 220 
km sec -1 with pf* = 0.34 GeV cm -3 and pg lax = 0.71 
GeV cm -3 ; iii) v Q = 270 km sec -1 with pf™ = 0.62 GeV 
cm -3 and p|f ax = 1.07 GeV cm -3 . 



B. Local fractional density 

The WIMP fractional density £ = pw/po is taken to 
be £ = min{l,r2 x /i 2 /(r2c , DM^ 2 )min}, in order to have 
rescaling [2^], when Q x h 2 turns out to be less than 
(flcDMh 2 )min (here (Qcdm h 2 )min is set to the value 
0.098, as derived at a 2er level from the results of Ref. 
0.) 



WIMP— nucleon cross section: hadronic 
uncertainties 



The ncutralino nucleon scalar cross section a. 



(nucleoli) 
scalar 



is mainly due to exchanges of the two CP-even neutral 
Higgs bosons, h and H, in the t-channel [26j and to 
squark-exchanges in the s- and u-channels [271 ] . The 
expression of cr^iax" 1 ^ may be found in [2|| . 

In the cross section representing the Higgs-exchanges, 
the factors involving the couplings between the Higgs 
bosons and the nucleon may be written as Ih.H = 
J2 q kq' H m q (N\qq\N) , where the coefficients k^' H depend 
on supersymmctric parameters, (N\qq\N) denotes the 
scalar density of the quark q in the nucleon, and m a is 
the quark mass. The coefficients k^' H are given in 



29] 
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The calculation of the quantities m q (N\qq\N) is usu- 
ally carried out by first expressing these quantities in 
terms of the pion-nucleon sigma term 

(JttN = x( m « + m d ) < N\uu + dd\N >, (4) 

of the quantity <7o, related to the size of the SU(3) sym- 
metry breaking, 

(T = -jijnu + m d ) < N\uu + dd - 2ss\N >, (5) 

and of the ratio r = 2m s /(m u + md). 

In fact, by assuming isospin invariance for quarks u 
and d, the quantities m q {N\qq\N) for light quarks may 
be written as 



41 MeV < ovat < 57 MeV 



(12) 



as derived from the pion-nucleon scattering amplitude, 
calculated at the so-called Cheng-Dashen point by Koch 
[39I ] . and from the evolution of the nucleon scalar form 
factor, as a function of the momentum transfer from t = 
2m\ to t = 0; evaluated in Ref. [40 1. 

In Ref. [29| we reconsidered the calculation of the co- 
efficients g u , gd in light of a new determination of aviv 
presented in Ref. [4lJ. In fact, the George Washing- 
ton University / TRIUM F g roup , using up-dated pion- 
nucleon scattering data [42j and a new partial-wave and 
dispersion relation analysis program, derived a range for 
0Vat [41 1 



m u < N\uu\N > ~ m d < N\dd\N >~ -a wN (6) 



m s < N\ss\N > ~ -r(a nN - cr ). 



(7) 



As for the heavy quarks c, b, t, one conveniently em- 
ploys the heavy quark expansion |30j to obtain 



m c < N\cc\N >~ m b < N\bb\N >~ m t < N\tt\N > 



2 

27 



m N - a nN + -r(a-xN - cr ) 



,(8) 



where mjv is the nucleon mass. 

In this way the quantities I^.h can be re-expressed as 



T _ h h - H _L L h > H 

lh,H — K u -tvpc9u + ^ 



(9) 



where 



g u ~ mi < N\ll\N > + 2 m h < N\hh\N > 
4 19 1 



ft/ 



to/ < 7V|ZZ|7V > + to s < iV|ss|iV > + 
- m h < N\hh\N > 



2 23 25 

^ ^(mAr + — ct^at + — r(cr 7rA r - cr )); (H) 

I stands for light quarks (I = u, d) and h denotes the 
heavy ones (h = c, b, t). 

The major problem here is that the three quantities 
&irN, 0o and r are all affected by sizable uncertainties 
which in turn dramatically affect the determination of 
the coefficients g u and gd- This problem was stressed in 
Ref. (29I . 31 1 (see also Refs. 32, 33 1 for earlier discussions 
of this point) and subsequently also recognized by other 
authors [H, Hi, M 



The range of a^N we considered in Ref. 3l| was 



55 MeV < a-xN < 73 MeV, 



(13) 



which turned out to be sizeably larger than the one of 
Eq. (|12p . Values of the nucleon scalar form factor at the 
Cheng-Dashen point higher than those of Ref. [i^] were 
also reported in Ref. [431 ] . 

In Ref. [2^] we stressed the dramatic importance that 
the uncertainties in the hadronic quantities ovat, cto and 
r have in a number of fundamental issues, namely: i) 
evaluation of WIMP detection rates in direct and also 
in some indirect searches (neutrino fluxes due to WIMP 
self-annihilation in Earth and Sun), ii) actual regions of 
the supersymmetric parameter space involved in searches 
for WIMPs, iii) connection between experimental event 
rates and relic abundance. 

It is unfortunate that the dichotomy between the two 
determinations in Eq. (TT2]) and in Eq. (TT3")) still persists. 
This fact is not only related to the experimental determi- 
nation of the relevant quantities in the pion-nucleon scat- 
tering but also to the intricacies involved in the derivation 
of ova/ from the experimental data. It is worth noting 
that the uncertainty inherent in the value of ovat could 
be even larger than the one exemplified by Eqs. (|12[) — (|T5j) 
(see Ref. [44J for further determinations). 

For definitcness, in the present paper we consider vari- 
ations of avAr in the range which is the union of the two 
intervals of Eqs. ([I^ ]) -([T5 ]) . i.e. 



41 MeV < Okn < 73 MeV. 



(14) 



In the presentation of our results, we will also consider 
a reference point, representative of a value of a^N which 
is within the narrow overlap of the two ranges of Eqs. 



The quantity ctq is taken in the range 45] 
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(T = 30 4- 40 MeV 



(15) 



and to the mass ratio r = 2m s /(m u + rrid) the default 
value r = 25 is assigned. 

We recall that the fractional strange-quark content of 
the nucleon y 



< N\ss\N > 



< N\uu + dd\N > ' 
is linked to ctq and a^N by the expression 



(16) 



y 



i - 



go 

O-kN 



(17) 



The reference point mentioned above, meant to repre- 
sent an estimate of the hadronic quantities in agreement 
with both ranges of Eq. (fT2")) and Eq. (fT3"|) . is defined by 
the values g u ^ e f = 123 MeV, gd,ref = 290 MeV (this set 
of values is the one employed also in our previous paper 
of Ref. 55j). As mentioned above, together with this 



representative point, we will also explicitly discuss the 
implications of the full range of a^jy given in Eq. lfT^ 



IV. THE ANNUAL MODULATION EFFECT 
MEASURED BY THE DAMA COLLABORATION 

The first results of the DAM A/LIBRA experiment on 
direct detection of dark matter particles have recently 
been presented Q. These data concern an exposure of 
0.53 ton yr. When added to the previous exposure of 
0.29 ton yr of the DAMA/Nal experiment 0], the to- 
tal exposure collected by the DAMA Collaboration, with 
the DAM A/LIBRA and the DAMA/Nal experiments to- 
gether, amounts to 0.82 ton yr. The analysis of the to- 
tal set of data shows an annual modulation effect in the 
event rate at 8.2 a C.L. This yearly modulation satis- 
fies all the features expected for an annual variation due 
to relic particles in our galactic halo [3] and is not ex- 
plained by systematic effects and/or seasonal variations 
of other various origins Q . 

The relic particles responsible for the DAMA annual 
modulation effect can be of many different kinds; a num- 
ber of possible candidates are recollected in Ref. (see 



e.g. 



also 
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49j, l50(). In the present paper we an- 
alyze the interpretation of the DAMA effect in terms of 
relic neutralinos. 

Indeed, we already invoked neutralinos 51| since the 
very first DAMA results of Ref. [H^] which showed an 
yearly variation in the direct signal. Our interpretation 
was further pursued, most notably in Refs. [9j, l53l . |54|. 




FIG. 1: t CT ( nuclc ° n) as a function of the WIMP mass. The 

3 SCcilciir 

region covered by a (red) slant hatching denotes the DAMA 
annual modulation region, under the hypothesis that the ef- 
fect is due to a WIMP with a coherent interaction with nuclei 
and without including the channeling effect. This region rep- 
resents the domain where the likelihood-function values differ 
more than 6.5 a from the null hypothesis (absence of modu- 
lation). It has been derived by the DAMA Collaboration by 
varying the WIMP galactic distribution function over the set 
considered in Ref. [19f | and by taking into account other uncer- 
tainties of different origins [57[ ■ The scatter plot represents su- 
persymmetric configurations calculated with the model sum- 
marized in Sect. [II] at the fixed representative set of values for 
the hadronic quantities characterized by: g u , re f = 123 MeV, 
gd,ref = 290 MeV. The (red) crosses denote configurations 
with a neutralino relic abundance which matches the WMAP 
cold dark matter amount (0.098 < Q x h 2 < 0.122), while the 
(blue) dots refer to configurations where the neutralino is sub- 
dominant (Q x h 2 < 0.098). The (blue) uniformly-shaded re- 
gion represents the extension of the scatter plot upwards and 
downwards, when the hadronic uncertainties reported in Eq. 
(Tl4)l are included (see text). 



where the focus was set on light neutralinos (i.e. neu- 
tralino with masses m x < 50 GeV). 

More recently, we presented a new investigation of the 
DAMA/Nal results [Hi], when the DAMA Collaboration 
discussed the possible implications of including the chan- 
neling effect in the analysis of their data 

In this paper we update our interpretation of the an- 
nual modulation effect in light of the experimental results 
improved by the quite sizable increase in total exposure 
of the recent DAMA experiments. 

From now on, when mentioning DAMA data we will 
mean the total set of data including both the DAMA/Nal 
experiment and the DAMA/LIBRA one, reported in 
Refs. S. 
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FIG. 2: £ a y™^on, ag a function of the WIMP mass The 
region covered by a (red) slant hatching denotes the DAMA 
annual modulation region, under the hypothesis that the ef- 
fect is due to a WIMP with a coherent interaction with nuclei 
and including the channeling effect. All other prerequisites of 
this region are as in Fig. [1] The scatter plot and the (blue) 
uniformly-shaded region are as in Fig. [T] 



In the following section we will discuss the annual mod- 



ulation regions, in the plane m x -^a. 



(nucleoli) 
calar 



denotes 



the WIMP mass), that have been derived by the DAMA 
Collaboration for the case of a WIMP with a coherent 
interaction with nuclei [j^ ] . The derivation has been car- 
ried out both in the case where the channeling effect is 
taken into account and in the one where channeling is not 
included. We recall that the features of the annual mod- 
ulation region in the plane m x - £ CT scaiar C ° n ^ depend sen- 
sitively on the specific model-dependent procedure em- 
ployed in the evaluation of the channeling effect. The 
regions reported hereby for the case of channeling cor- 
respond to the case where channeling is included with 
the model explained in Ref. Q. Actually, the extent 
by which the channeling effect occurs when a putative 
WIMP traverses a Nal crystal is still under study. For 
this reason, in our analysis we consider both cases of no- 
channeling and of channeling with the model of Ref. [?J • 
One expects that the actual physical situation is com- 
prised within these two cases. 



A. Annual modulation regions (convolution over a 
class of distribution functions) 

We start our analysis by showing in Figs. QJ[2] the 
DAMA annual modulation regions in a plot of ^c^ilr 
versus the WIMP mass. These have been derived by the 



DAMA Collaboration for the case of a WIMP with a co- 
herent interaction with nuclei [56[ , by varying the WIMP 
galactic distribution function (DF) over the set consid- 
ered in Ref. [191 ] and by taking into account other uncer- 
tainties of different origins [571 ] ■ The DAMA regions are 
denoted by a (red) slant hatching; they represent regions 
where the likelihood-function values differ more than 6.5 
a from the null hypothesis (absence of modulation). Fig. 
[T] refers to the case in which the channeling effect is not 
included, whereas Fig. [2] displays the case where chan- 
neling is included. 

Figs. [TH2] also show the (blue) uniformly-shadowed 
region, which represents the physical ncutralino region 
as derived within our effective MSSM. The scatter plot, 
common to Fig.l and Fig. 2, denotes the results of our 
evaluations, when a scanning over the supcrsymmetric 
parameter space is performed, at the fixed representa- 
tive set of values for the hadronic quantities mentioned 
in Sect. IIII CI This representative set is characterized by 
the values: g u ,ref = 123 MeV, gd.ref = 290 MeV (we re- 
call that this set is the one employed also in our previous 
paper of Ref. (EsJ). The (red) crosses and the (blue) dots 
of the scatter plot denote configurations with no-rescaling 
and those with rescaling of the local density, respectively 
(see Sect. IIII Bp . The features of the supcrsymmetric 
configurations belonging to the scatter plot are discussed 
in detailed in Ref. [9||. 

The uniformly-shaded region displayed in Figs. HH2] 
represents the extension of the scatter plot upwards and 
downwards, when the hadronic uncertainties reported in 
Eq. (fl4j) are included. The range of the pion-nucleon 
sigma term of Eq. (|13| is responsible for the upper ex- 
tension of the physical region, as compared to the repre- 
sentative scatter plot, by an enhancement factor of about 
2-3, whereas the range of Eq. lfl"2|) generates the lower 
extension by a suppression factor of order 8-9. These 
numbers follow immediately from the formulae in Sect. 
IIII CI by taking into account that the dominant term 
in the quantity Ih,h is the one involving gd- Thus, the 
scatter plot for any value set of hadronic quantities with 
given values of gd and g u can approximately be obtained 
from the one corresponding to the reference set of val- 
ues, characterized by g u ,ref = 123 MeV, gd,ref = 290 
MeV, by scaling the reference scatter plot by the fac- 
tor (gd/ gd,ref) 2 ■ However, notice that, in deriving the 
boundaries of the full theoretical region in Figs. HH21 the 
full expression of Eq. © has been used. 

From Figs. 1-2 it is clear that the DAMA annual mod- 
ulation region is largely compatible with the theoretical 
predictions for relic neutralinos with masses m x < 100 
GeV, in particular for neutralinos within the low-energy 
funnel for m x < 50 GeV. This occurs, whether or not the 
channeling effect is included. 
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B. Annual modulation regions for single halo 
models 



We turn now to the analysis of the annual modulation 
regions for specific forms of the WIMP distribution func- 
tion. First we discuss in detail our reference model, the 
cored isothermal sphere, mentioned in Sect. IIII Al (de- 
noted as Evans logarithmic model, or Al model, in Ref. 
[lit); we will comment about some other DFs afterwords. 

FigsiSHU display the theoretical predictions of our su- 
persymmetric model (already shown in FigsHHSJ) to- 
gether with the DAMA annual modulation regions [j^ . 
under the hypothesis that the WIMP-nucleus interaction 
is coherent and that the velocity DF is given by the cored 
isothermal sphere. The various insets refer to the repre- 
sentative values of vq and po discussed before, that is: 
i) v = 170 km sec" 1 with pf lu = 0.20 GeV cm" 3 and 
p max = q 42 GeV cm" 3 ; ii) v = 220 km sec" 1 with 
fi? in = 0.34 GeV cm" 3 and p^ nax = 0.71 GeV cm" 3 ; iii) 
V = 270 km sec" 1 with p^ in = 0.62 GeV cm" 3 and 
p™ = 1.07 GeV cm -3 . v esc is set to the value v esc = 
650 km sec -1 . 

These DAMA annual modulation regions represent do- 
mains where the likelihood-function values differ more 
than 6.5 (J from the null hypothesis (absence of modula- 
tion). The notations for the various regions and for the 
scatter plot are the same as in Figs. [TH2J Fig. prefers to 
the case in which the channeling effect is not included, 
and Fig. |4]to the case where channeling is included. 

It is remarkable that relic neutralinos are able to pro- 
vide a good fit to the experimental data. In case of no- 
channeling, low values of vq and po (^o — 170 km sec -1 
and po — 0-2 GeV cm -3 ) appear to be somewhat disfa- 
vored, though in this case neutralinos with m x ~ 60-100 
GeV could be involved. The agreement between exper- 
imental data and theoretical evaluations in our model 
looks very good for the other combinations of vq and po 
values, with an overall preference for neutralinos of low 
mass. For the case where the channeling is included ac- 
cording to the modelling of Ref. Q, experimental data 
favor values of vq and po which are in the low-medium 
side of their own physical ranges, i.e. v ~ (170 - 220) 
km sec -1 and po ~ (0.3 — 0.4) GeV cm -3 and neutralino 
masses in the mass range m x ~ (7 — 30) GeV. 

When no rotation of the halo is considered, the features 
of the annual modulation region in the m x — Ccr"^ 1 ™" 
plane do not differ much when the galactic DF is varied 
among many of the galactic DFs considered in Ref. [3] . 
For instance, for a matter density with a Navarro-Frcnk- 
Whitc profile (A5 model of Ref. [19() or for an isothermal 
model with a non- isotropic velocity dispersion (Bl model 
of Ref. [lji]) the physical situations are very similar to 
the ones depicted in Fig. 3-4. 



However, in the case of DFs with triaxial spatial distri- 
butions (within the class D of Ref. [lj^) and for models 
with a co-rotating halo there can be an elongation of the 
annual modulation region towards heavier masses (these 
are generic characteristic features which can be derived 
from the analysis of Ref. 19(). These features are dis- 



played in Figs. 1-2 for masses up to the value of 100 GeV, 
which corresponds to the upper extreme of the annual- 
modulation region provided up to now by the DAMA 
Collaboration. 

We wish to recall that the distribution of WIMPs in the 
Galaxy could deviate from the models mentioned above, 
mainly because of the presence of streams. For modifica- 
tion of the annual modulation region in these instances 
see Ref. 0, [58 1 . 



V. UPPER BOUNDS ON THE 
WIMP-NUCLEON SCALAR CROSS SECTION 

Other experiments of WIMP direct detection, different 
from DAMA/Nal and DAMA/LIBRA, do not currently 
have the capability of measuring the annual modulation 
effect; they can only provide upper bounds for the ex- 
pected signals 0, 0, S @| • O nce a specific form for the 
WIMP DF is taken and the relevant parameters fixed, 
these bounds can be converted into constraints on the 
WIMP-nucleus cross section. It is obvious that, in order 
to derive solid constraints on the WIMP physical prop- 
erties, a conservative approach has to be followed in the 
selection of the WIMP distribution function; i.e. among 
the wide variety of DFs, one has to select the ones which 
imply the smallest responses in the detector. 

As an example, we take the data recently published 
by the CDMS Collaboration From these data, ap- 
plied to a WIMP with a coherent interaction with nu- 
clei, one can derive the upper bounds for the quantity 
£ CT scaiar C ° n ' displayed in Fig. [5j for a number of vari- 
ous DFs. The strong dependence of the upper bounds 
on the assumed DF is apparent. Superimposed to these 
limits is our theoretical region for relic neutralinos. We 
note that the conservative upper bound, established by 
the Bl distribution function, does not exclude any of the 
light neutralino configurations, when the hadronic uncer- 
tainties are taken into account (only some configurations 
with neutralino masses above 60 GeV are excluded). Of 
course, other WIMP distribution functions, such as the 
isothermal sphere, would introduce constraints on the su- 
persymmetric population. 

However, one has to notice that the aforementioned 
bounds are obtained through involved procedures for 
discriminating electromagnetic signals from recoil events 
and through delicate subtractions meant to separate pu- 
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FIG. 3: C a scaiar° 0n ' as a function of the WIMP mass. The region covered by a (red) slant hatching denotes the DAMA annual 
modulation region, under the hypothesis that the effect is due to a WIMP with a coherent interaction with nuclei and without 
including the channeling effect. This region represents the domain where the likelihood-function values differ more than 6.5 a 
from the null hypothesis (absence of modulation) . It has been derived by the DAMA Collaboration by assuming that the WIMP 
distribution function is given by the cored isothermal sphere (denoted as Evans logarithmic model, or Al model, in Ref. 
and using the parameters of set A of Sect. 7.2 of the first paper of Ref. (jj. The scatter plot and the (blue) uniformly-shaded 
region are as in Fig. [T] 



tative WIMP signals from neutron-induced events. A 
major critical point in these experiments and related 
analyses is that the very signature (the annual modu- 
lation) of the searched signal cannot be employed in ex- 
tracting the authentic events. Possible problems of sta- 
bility in the acquisition parameters can affect the rejec- 
tion procedures applied to a large number of events, as 
well as the determination of the threshold and of the en- 



ergy scale [59(. In view of possible sizable uncertainties 



involved in these procedures, we conservatively do not 
implement the upper bounds discussed in this section, 
while comparing theoretical expectations for relic neu- 
tralinos to the annual modulation data of Ref. [3] ■ 

It is also worth noting that the upper bounds of Ref. 
0, even when taken at their face value, are however not 
in conflict with the annual modulation data and the theo- 
retical neutralino interpretation for masses ^7—10 GeV. 
Compatibility of the DAMA/LIBRA annual-modulation 
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FIG. 4: C a scaiar° 0n ' as a function of the WIMP mass. The region covered by a (red) slant hatching denotes the DAMA annual 
modulation region, under the hypothesis that the effect is due to a WIMP with a coherent interaction with nuclei and including 
the channeling effect. All other prerequisites of this region are as in Fig. [3] The scatter plot and the (blue) uniformly-shaded 
region are as in Fig. [1] 



results with the upper bounds of other experiments for 
coherently-interacting WIMPs with the masses < 10 
GeV is also found in two papers 6fJ, IfUj , appeared con- 



VI. GALACTIC ANTIPROTONS 



comitantly with the present article. 



Finally, we recall that another experiment (KIMS [6(]), 
running at present with a detector of about 104 kg of 
Csl crystals, is meant to provide a measurement of the 
WIMP annual modulation in the future. 



In Sects. IIV AlllVBl we have seen that the agreement 
of the DAMA data with theoretical evaluations for relic 
neutralinos is quite good for light ncutralino masses and 
for wide intervals of the astrophysical quantities po and 
vq within the ranges discussed in Sect. IIII Al However, 
it turns out that the data of cosmic antiprotons fluxes 
can put stringent limits on light ncutralino configurations 
55|, Ell- 



in fact, in Ref. |63[ it is shown that the experimental 
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FIG. 5: The solid lines show the upper limit on the quantity 
t (micleon) ag a function of the WIM p mass m v for the CDMS 

~ scalar a. 

detector [5J and for v CBC = 650 km sec -1 . The (red) median 
line refers to the standard isothermal sphere with vq = 220 
km sec -1 and p = 0.3 GeV cm -3 (model AO of Ref. 0). 
The (black) upper and lower curves refer to model Bl with 
vo = 170 km sec -1 (upper solid line) and model C3 with vo = 
270 km sec -1 (lower solid line). The short-dashed line refers 
to model C3 with maximal counter-rotation of the galactic 
halo. The long-dashed lines show the upper limits for CDMS 
in the case of a lower escape velocity v eBC = 450 km sec -1 : 
the upper line refers to model Al, the lower one to model 
C3. For model Bl, the limit coincides with the corresponding 
solid line. The scatter plot and the (blue) uniformly-shaded 
region are as in Fig. [1] Other specifications in the text. 



antiproton spectrum is fitted very well by the secondary 
component from cosmic rays spallation, calculated with 
the set of the diffusion parameters which is derived from 
the analysis of the boron-to-carbon ratio (B/C) compo- 
nent of cosmic rays 



64 



Indeed the calculated secondary 
flux of cosmic antiprotons fits the experimental data with 
with a x 2 = 33.6 with 32 data points, and with an un- 
certainty of 20%. This means that very little room is left 
for a possible primary component of antiprotons gener- 
ated by an exotic origin (neutralino self-annihilation in 
our case). 

Let us recall that in Ref. the secondary antiproton 
spectrum, generated by spallation processes, was prop- 
agated using a two-zone diffusion model described in 
terms of five parameters. Two of these parameters, Kq 
and 5, enter the expression of the diffusion coefficient: 
K = Kof3R 5 (R is the particle rigidity); the other three 
parameters are the Alfen velocity Va, the velocity of the 
convective wind V c , and the thickness L of the two large 
diffusion layers which sandwich the thin galactic disk 64] . 



case 


S 


Ko 


L 


V c 


Va 






[kpc 2 /Myr] 


[kpc] 


[km s -1 ] 


[km s -1 ] 


max 


0.46 


0.0765 


15 


5 


117.6 


med 


0.70 


0.0112 


4 


12 


52.9 


min 


0.85 


0.0016 


1 


13.5 


22.4 



TABLE I: Astrophysical parameters of the two-zone diffusion 
model for galactic cosmic rays propagation, compatible with 
B/C analysis and yielding the maximal, median and minimal 
primary antiproton flux [651 ] - 



When studying the primary antiprotons one can extract 
the following three sets of propagation parameters: the 
best-fit (on B/C) set (denoted as median), together with 
the sets which yield the minimal and the maximal pri- 
mary antiproton fluxes [6^]. The values of these three 
sets are given in Table I. 

We proceed now to analyze the extent of compatibil- 
ity of the neutralino configurations which fit the DAMA 
results with the present data on cosmic antiprotons. 
Among the six sets of values for the parameters vq and po 
analyzed in Figs. [3H1J let us consider the following ones 
(the same already considered in Ref. 0): A) i> = 170 
km sec -1 , pg lin = 0.20 GeV cm -3 ; B) v = 170 km 
sec -1 , pg lax = 0.42 GeV cm -3 ; C) v = 220 km sec -1 , 

„min 

Po 



= 0.34 GeV cm -3 . 
In Figs. EH?] we give the antiproton fluxes at p kinetic 
energy Tp = 0.23 GeV, as a function of the neutralino 
mass for a cored isothermal halo and for the neutralino 
configurations selected by the DAMA regions shown in 
Fig. [3HU respectively. Fig. prefers to the case in which 
channeling is not included in the derivation of the DAMA 
regions, Fig. [7] to the case with channeling included. 

At variance with what displayed in Figs. HHH where 
the scatter plot was evaluated at the reference point: 
g u , ref = 123 McV, g d<ref = 290 MeV, in Figs. EHT] the 
scatter plots are the results not only of the scan over the 
supersymmctric parameter space but also of the varia- 
tions of the quantities g u and gd, as given by Eqs. (|10p - 
(jlip . when ovtv and do are varied in the ranges of Eq. 
(fl"4|) and Eq. (|15|) . respectively, and r is put at the default 
value r = 25. 

From the results shown in Figs. [5H7] we see that, 
though a number of configurations arc excluded by the 
BESS data 66j, many others are perfectly compatible 



with BESS and in principle accessible to PAMELA [67 1 
and AMS-02 [gl]. More specifically, for set A, most of 
the neutralino configurations are unconstrained by the 
galactic antiproton data, except for a group of them in 
the case of the maximal set of the diffusion parameters 
and when channeling is included; a sizable number of 
configurations are at the level of possible investigation. 
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FIG. 6: Antiproton flux, at p kinetic energy Tp = 0.23 GeV, generated by the neutralino configurations selected by the DAMA 
data, when a cored isothermal halo is employed and the channeling effect is not included. Each row corresponds to a different 
set of cosmic rays propagation parameters: the upper, median and lower rows refer to the sets which provide the maximal, 
median and minimal antiproton flux, respectively (see Table I). The three columns refer to the sets of vo and po values denoted 
as sets A, B and C in the text. The bold (colored) points refer to configurations compatible with the DAMA regions for 
the case in which the channeling effect is not included (see Fig. [3]). In selecting the allowed configurations the hadronic 
uncertainties have been taken into account. The points are differentiated as follows: (red) crosses denote configurations with 
a neutralino relic abundance which matches the WMAP cold dark matter amount (0.098 < Q. x h 2 < 0.122), while (blue) dots 
refer to configurations where the neutralino is subdominant (Jl x h 2 < 0.098). The light gray points denote configurations with 
a neutralino-nucleon scattering cross section outside the corresponding DAMA allowed region. The solid horizontal line shows 
the maximal allowable amount of antiprotons in the BESS data [g^] over the secondary component; the dashed and dotted lines 
denote estimates of the PAMELA [63] and AMS-02 [73|] sensitivities to exotic antiprotons for 3 years missions, respectively. 



Sets B and C, due to their corresponding higher values of 
Po, are more sensitive to the p constraints (but also, to a 
large extent, accessible to PAMELA and AMS), though 
prevalently for sets of the diffusion parameters close to 



the maximal set. 

It is worth noticing explicitly that for the other sets 
of v and p discussed in Sect. IIVBI and in Figs. [3HH 
but not considered here, one would obtain plots similar 
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FIG. 7: Antiproton flux, at p kinetic energy Tp = 0.23 GeV, generated by the neutralino configurations selected by the DAMA 
data, when a cored isothermal halo is employed and the channeling effect is included. Notations are as in Fig. [6] except that 
here the neutralino configurations of the scatter plot are those selected by the DAMA regions of Fig. [4] 



to the ones displayed in Figs. [6H3 with scatter plots 
rescalcd according to the power p§; thus for these sets the 
p constraints would be more severe than in the previous 
cases. 



powerful in providing stringent constraints, are some- 
what problematic in providing positive signals for exotic 
production. Different is the case for antidcuterons which 
we discuss in the following section. 



As we noticed above, a sizable number of neutralino 
configurations are at the level of the sensitivities of cur- 
rent experiments on cosmic antimatter. However, for 
the reasons explained above, the corresponding primary 
fluxes would be rather difficult to be disentangled from 
the secondary flux (notice that primary and secondary 
fluxes have also a verysimilar behavior as functions of 
the p kinetic energy) 65[. The cosmic antiproton data, 



VII. SEARCH FOR ANTIDEUTERONS IN THE 
GALACTIC HALO 

In Ref. [6^ it was shown that the antideuteron spec- 
tra derived from DM self-annihilation are much flatter 
than the expected astrophysical component for kinetic 
energies Tq < 2-3 GcV/n. Recently, the calculation of 
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FIG. 8: Antideuteron flux, at D kinetic energy Tg = 0.23 GeV/n, generated by neutralino configurations selected by the 
DAMA data when a cored isothermal halo is employed and the channeling effect is not included. Each row corresponds to 
a different set of cosmic rays propagation parameters: the upper, median and lower rows refer to the sets which provide the 
maximal, median and minimal antiproton flux, respectively (see Table I). The three columns refer to the sets of vo and po values 
denoted as sets A, B and C in the text. The bold (colored) points refer to configurations compatible with the DAMA regions 
for the case in which the channeling effect is not included (see Fig. [3]). In selecting the allowed configurations the hadronic 
uncertainties have been taken into account. The points are differentiated as follows: (red) crosses denote configurations with 
a neutralino relic abundance which matches the WMAP cold dark matter amount (0.098 < £l x h 2 < 0.122), while (blue) 
dots refer to configurations where the neutralino is subdominant (Q x h 2 < 0.098) The light grey (green on-line) points denote 
supersymmetric configurations yielding an exceedingly high antiproton flux (see Fig. |6}. The horizontal lines refer to estimated 
sensitivities to antideuterons of the GAPS (dashed) and AMS (dotted) detectors. 



the primary and secondary antideuteron fluxes has been 
performed [7o| in the framework of the full propagation 
model outlined in the previous section, encoding all the 
possible uncertainty sources. InRef. 0, |7l| it has been 
shown that antideuterons offer a very promising signal for 



the indirect detection of intermediate and low mass DM 
particles by means of future detectors such as GAPS on 
long and ultra-long duration balloon (LDB and ULDB) 
flights 0] and AMS-02 for three years data taking 

In the present Section, we present results on the fluxes 
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FIG. 9: Antideuteron flux, at D kinetic energy Tg = 0.23 GeV/n, generated by neutralino configurations selected by the 
DAMA data when a cored isothermal halo is employed and the channeling effect is included (see Fig. [4| . The light grey (green 
on-line) points denote supersymmetric configurations yielding an exceedingly high antiproton flux (see Fig. [7]l . Other notations 
are as in Fig. [8] except that here the neutralino configurations of the scatter plot are those selected by the DAMA regions of 
Fig. H 



of antideuterons produced from DM self-annihilation in 
the galactic halo following the procedure explained in 
Ref. [n3] (we refer to this paper for all details). In Fig. 
8 we display the antideuteron flux calculated at D ki- 
netic energy of Tq = 0.23 GeV/n as a function of the 
neutralino mass for the cored isothermal halo Al of Eq. 
J3]) . The scatter plot represents the supersymmetric con- 
figurations selected by DAMA data when the channel- 
ing effect is not included, while Fig. 9 shows the same 
quantities when the channeling effect is included. The 
configurations have been selected with the same method 



as in the previous Section. The three upper, median and 
lower figures correspond to the maximum, median and 
minimum set of propagation parameters of Table I. The 
dotted and dashed horizontal lines assess the sensitivities 
of AMS-02 and GAPS ULDB, respectively. The light 
gray (green) points correspond to supersymmetric con- 
figurations yielding an exceedingly high antiproton flux 
(see Figs. [BH7]). Figs. [SHE show that a sizable number 
of neutralino configurations compatible with the annual 
modulation data can generate signals accessible to an- 
tideuteron searches planned for the next years. It is no- 
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FIG. 10: Solid line: antideuteron flux for a DM halo of m x = 
20 GeV neutralinos compatible with the DAMA modulation 
effect. Dashed line: secondary component [7(J. Propagation 
parameters are set at the median configuration of Table I. The 
upper dashed horizontal line is the present BESS upper limit 
on the search for cosmic antideuterons. The three horizontal 
solid (blue) lines are the estimated sensitivities for (from top 
to bottom): AMS-02, GAPS on a LDB and GAPS on an 
ULDB. 



tabic that a bunch of configurations not accessible to an- 
tiproton signals could instead be probed by antideuteron 
searches. 

In Fig. 1101 by way of example, we show the an- 
tideuteron flux for a DM halo made of m x = 20 GeV 
neutralinos compatible with the DAMA effect, with po = 
0.34 GeV cm -3 and setting the propagation parameters 
to the median values. The lower, dashed line corresponds 
to the secondary antideuteron flux as obtained in Ref. 
[ToT ] . Fluxes are modulated at solar minimum. The three 
horizontal lines arc the estimated sensitivities for (from 
top to bottom): AMS-02, GAPS LDB and GAPS ULDB 
flights. The primary flux stands well above the back- 
ground and exceeds the experimental reach of next gen- 
eration detectors. The result is noticeable: both GAPS 
and AMS-02 will have the capability of clearly detecting 
antideuterons produced from a neutralino halo compati- 
ble with the positive signal found in DAMA data. 



VIII. MEASUREMENTS AT LHC 

Naturally, the viability of the interpretation of the 
DAMA annual modulation data in terms of relic neu- 
tralinos depends ultimately on the results which LHC 
will provide on supersymmetric theories. A thorough in- 
vestigation on the search for light neutralinos at LHC has 
been recently carried out in Ref. 74] . In this paper two 



specific scenarios are analysed, both of which dictated by 
cosmological properties pertaining to light neutralinos. 

A first scenario, scenario A, is identified by the follow- 
ing sector of the supersymmetric parameter space: Mi ~ 
10 GeV, \fj,\ ~ (100- 200) GeV, m A ~ 90 GeV, tan/3 - 30 
-45,-1<j4<+1; the other supersymmetric parameters 
are not a priori fixed. In this scenario the cosmologi- 
cal bound il x h 2 < (ricr>M^ 2 )max is satisfied because of 
a neutralino self annihilation mediated by the light A— 
boson, combined with a high value of tan j3 and a sizeable 
bino-higgsino mixing in the neutralino composition. 

A second scenario, Scenario B, is identified by the fol- 
lowing sector of the supersymmetric parameter space: 
Ml - 25 GeV, \fi\ > 500 GeV, tan/3 < 20; mj > (100 
- 200) GeV, -2.5 < A < +2.5; the other supersymmet- 
ric parameters are not a priori fixed. In this case the 
cosmological constraint il x h 2 < (^lcDMh 2 ) max is satis- 
fied because of stau-exchange contributions (in the t, u 
channels) to neutralino self-annihilation cross section. 

In Ref. 74| the signals expected at LHC in the two 



scenarios are discussed through the main (sequential and 
branched) chain processes, started by a squark produced 
in the initial proton-proton scattering. Branching ratios 
and the expected total number of events are derived in 
the various benchmarks defined within the two scenarios. 



On the basis of these results it is proved in Ref. 74 1 
that LHC should allow an efficient exploration of the su- 
persymmetric parameter regions compatible with light 
neutralinos. Due to the characteristic features of these re- 
gions, the measurements of LHC should be able to prove 
or disprove the supersymmetric model considered in the 
present paper and then validate or not our interpretation 
of the annual modulation effect in terms of relic neutrali- 



IX. CONCLUSIONS 

In this paper we have analyzed the most recent exper- 
imental data on direct searches for dark matter parti- 
cles in the galactic halo. We have discussed the various 
features involved in an interpretation of these data in 
terms of relic neutralinos; the supersymmetric scheme 
employed is an effective MSSM scheme at the elec- 
troweak scale, without gaugino-mass unification at a 
Grand Unified scale. The role of the uncertainties af- 
fecting the neutralino-quark couplings, because of the 
involved hadronic quantities, are critically discussed and 
included in our considerations. 

First we have considered in detail the results of the 
DAMA Collaboration which, by a combined analysis of 
the DAMA/Nal and the DAM A/LIBRA experiments, for 
a total exposure of 0.82 ton yr, provide evidence of an 
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annual modulation effect at 8.2 a C.L. Comparison of our 
theoretical evaluations with the DAMA data has been 
carried out both in the case in which channeling effect 
is included in the DAMA analysis and in the one where 
it is excluded. The DAMA results presented here refer 
uniquely to the case in which the putative WIMP has a 
coherent interaction with the nuclei in the detector and 
are therefore translated in terms of regions in the m x — 

£ CT scaiar COn ^ pl ane i f° r different selections of the galactic 
distribution functions [56|]. We wish also to recall that 
DM-detector interaction mechanisms or DM candidates 
different from the one considered in the present paper 
could be the origin of the annual modulation effect [2|. 

The roles the hadronic uncertainties and of the effect 
of channeling are kept separated at any stage of our dis- 
cussion on DAMA data; the analysis has been carried 
out in such a way that, once hopefully some of the afore 
mentioned uncertainties are resolved in the future, one 
can easily employ our present results to narrow down the 
physical relevant regions. 

By considering first the DAMA data when a convolu- 
tion of galactic distribution functions is considered, we 
have shown that the annual modulation region is largely 
compatible with the theoretical predictions for relic neu- 
tralinos with masses m x < 100 GeV, and in particular 
for neutralinos within the low-energy funnel for m x < 50 
GeV. This occurs, whether or not the channeling effect 
is included. 

We have then pursued our analysis by employing the 
cored isothermal sphere DF and shown that relic neu- 
tralinos fit quite well the annual modulation data. For 
the case where the channeling is included according to 
the modeling of Ref . Q , experimental data favor values 
of vq and po which are in the low-medium side of their 
own physical ranges, i.e. vq ~ (170 - 220) km sec -1 and 
Po — (0.3 — 0.4) GeV cm -3 and neutralino masses in the 
mass range m x ~ (7 — 30) GeV. In case of no-channeling, 
low values of ^o and po (vq ~ 170 km sec -1 and po ~ 0.2 
GeV cm -3 ) appear to be somewhat disfavored, though 
in this case neutralinos with m x ~ 60-100 GeV could be 
involved. The agreement between experimental data and 
theoretical evaluations looks very good for the other com- 
binations of Vq and po values, with an overall preference 
for neutralinos of low mass. We have also commented 
about the results one would obtain when some other DF 
is selected. 

The neutralino populations selected on the basis of the 
annual modulation data have been analyzed in terms of 
the antiproton fluxes which they would produce in our 
halo. We have shown that many of them arc fully com- 
patible with the current stringent bounds on cosmic an- 
tiprotons, especially for values of local dark matter den- 
sity po and local rotational velocity vq in the low side 



of their physical ranges, and for values of the diffusion 
parameters not too close to the values of their maximal 
set. 

It was also derived that forthcoming measurements of 
galactic antidcuterons will have good chances to explore 
fluxes generated by the self-annihilation of the neutralino 
configurations compatible with the annual modulation 
data. 

The upper bounds on C CT icaiar° n ' ) ' derived by other ex- 
periments of WIMP direct searches 0, 0, 0, @], have 
been analyzed by considering various WIMP galactic DFs 
and by including the uncertainties induced by the previ- 
ously mentioned hadronic quantities. It is concluded that 
no conservative upper bound can now exclude the light 
(m x < 50 GeV) neutralino population within the model 
discussed in Sect. [TTJ 

We have pointed out that LHC is expected to allow 
an efficient exploration of the supersymmetric parameter 
regions compatible with light neutralinos and the annual 
modulation data. Due to the characteristic features of 
the physical regions involved, the measurements of LHC 
could be able to prove or disprove the supersymmetric 
model considered in the present paper and then validate 
or not our interpretation of the annual modulation effect 
in terms of relic neutralinos. 

It is finally worth noting that the relic neutralinos of 
light masses, considered in the present paper and pre- 
viously in Refs. 0, Hi)], can be of relevance also in 
many interesting astrophysical contexts different from 
the ones discussed before. For instance, by their self- 
annihilations, these neutralinos may be the origin of 
electron-positron pairs which could interact with pho- 
tons of the cosmic microwave back grou nd producing a 
peculiar Sunyaev-Zeldovic effect 75|, 76|. Also, positrons 
and electrons created by annihilation of light neutralinos 
in the galactic center might generate a synchrotron radia- 
tion responsible for an excess of microwave emission from 
a region within ~ 20° of the galactic center [xtJ (see also 
Ref. ) . Light neutralino dark matter could also imply 



a sizable primordial 6 Li abundance 79j, a feature partic- 
ularly interesting in view of recent new determinations of 
this abundance in metal-poor halo stars 80 1 . 
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